Curvature-enhanced spin-orbit coupling in a carbon nanotube 
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Structure of the spin-orbit coupling varies from material to material and thus finding the correct 
spin-orbit coupling structure is an important step towards advanced spintronic applications. We 
show theoretically that the curvature in a carbon nanotube generates two types of the spin-orbit 
coupling, one of which was not recognized before. In addition to the topological phase-related 
contribution of the spin-orbit coupling, which appears in the off'-diagonal part of the effective Dirac 
Hamiltonian of carbon nanotubes, there is another contribution that appears in the diagonal part. 
The existence of the diagonal term can modify spin-orbit coupling efi'ects qualitatively, an example of 
which is the electron-hole asymmetric spin splitting observed recently, and generate four qualitatively 
different behavior of energy level dependence on parallel magnetic field. It is demonstrated that the 
diagonal term applies to a curved graphene as well. This result should be valuable for spintronic 
applications of graphitic materials. 
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I. INTRODUCTION 

Graphitic materials such as carbon nanotubes (CNTs) 
and graphenes are promising materials for spintronic ap- 
plications. Various types of spintronic devices arc re- 
ported such as CNT-based three terminal magnetic tun- 
nel junctionsi, spin diodes^, and graphene-based spin 
valves^. Graphitic materials are believed to be excellent 
spin conductors^. The hyperfine interaction of electron 
spins with nuclear spins is strongly suppressed since ^^C 
atoms do not carry nuclear spins. It is estimated that 
the spin relaxation time in a CNT— and a graphene^ is 
limited by the spin-orbit coupling (SOC). 

Carbon atoms are subject to the atomic SOC Hamil- 
tonian Hso- In an ideal flat graphene, the energy shift 
caused by Hgo is predicted to be ~ 10"'^ meV— i^. Re- 
cently it is predicted^i^ that the geometric curvature 
can enhance the effective strength of the SOC by orders 
of magnitude. This mechanism applies to a CNT and 
also to a graphene which, in many experimental situa- 
tions, exhibits nanometer-scale corrugations^^. There is 
also a suggestion' that artificial curved structures of a 
graphene may facilitate device applications. 

A recent experiments^ on ultra-clean CNTs measured 
directly the energy shifts caused by the SOC, which 
provides an ideal opportunity to test theories of the 
curvature-enhanced SOC in graphitic materials. The 
measured shifts are in order-of-magnitude agreement 
with the theoretical predictions^, confirming that the 
curvature indeed enhances the effective SOC strength. 
The experiment revealed discrepancies as well; While ex- 
isting theories predict the same strength of the SOC for 
electrons and holes, which is natural considering that 
both the conduction and valence bands originate from 
the same tt orbital, the experiment found considerable 
asymmetry in the SOC strength between electrons and 
holes. This electron- hole asymmetry implies that existing 
theories of the SOC in graphitic materials are incomplete. 

In this paper, we show theoretically that in addition 



to effective SOC in the off-diagonal part of the effective 
Dirac Hamiltonian, which was reported in the existing 
theories^i^, there exists an additional type of the SOC 
that appears in the diagonal part both in CNTs and 
curved graphenes. It is demonstrated that the combined 
action of the two types of the SOC produces the electron- 
hole asymmetry observed in the CNT experimenti^ and 
gives rise to four qualitatively different behavior of en- 
ergy level dependence on magnetic field parallel to the 
CNT axis. 

This paper is organized as follows. In Sec. |TT1 we show 
analytical expressions of two types of the effective SOC 
in a CNT and then explain how the electron-hole asym- 
metric spin splitting can be generated in semiconducting 
CNTs generically. Section Hill describes the second-order 
perturbation theory that is used to calculate the effective 
SOC, and tight-binding models of the atomic SOC and 
geometric curvature. Section HVl reports four distinct en- 
ergy level dependence on magnetic field parallel to the 
CNT axis. We conclude in Sec. |V] with implications of 
our theory on curved graphenes and a brief summary. 



II. EFFECTIVE SPIN-ORBIT COUPLING IN A 

CNT 

We begin our discussion by presenting the first main 
result for a CNT with the radius R and the chiral angle 
e {0<e<n/6, 0(7r/6) for zigzag (armchair) CNTs). We 
find that when the two sublattices A and B of the CNT 
are used as bases, the curvature-enhanced effective SOC 
Hamiltonian Tisoc near the K point with Bloch momen- 
tum K becomes 



{SyR)ay {5'jR)ay 



(1) 



where Uy represents the real spin Pauli matrix along the 
CNT axis. The pseudospin is defined to be up (down) 
when an electron is in the sublattice A (B). Here the off- 
diagonal term that can be described by a spin-dependent 
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topological phase are reported in Refsi^i^ but the diago- 
nal term was not recognized before. Expressions for the 
parameters Sk and 5'^^ are given hy^ 



Sk 
R 



A,oa(£s-£p)(V7p + F/p) 



R 



and 



cos W 



R 2V3iV--VJl) R 



(2) 



(3) 



where Ago ^ 12meV— is the atomic SOC constant, a is 
the lattice constant 2.49A, and es(p) is the atomic energy 

for the s{p) orbital. Here, Vfp and VpJ-'^'' represent the 
coupling strengths in the absence of the curvature for the 
a coupling between nearest neighbor s and p orbitals and 
the 7r((T) coupling between nearest neighbor p orbitals, 
respectively. Note that the \S'^\ has the ^-dependence, 
whose implication on the CNT energy spectrum is ad- 
dressed in Sec. |lVl For K' point with K' = -K, H^^^ is 



^soc 

given by Eq. ^ with Sk and 6'^ replaced by 6k' = —^k 
and S'^, = — (^K' respectively. 

Implications of the diagonal term of the SOC be- 
come evident when Eq. ([1]) is combined with the two- 
dimensional Dirac Hamiltonian Tioirac of the CNT. For 
a state near the K point with the Bloch momentum K + k 
[k = ik^,ky), |k| < |K|], Hoirac becomesi^ 



CNT, the diagonalization of 'H^^j.^^+Ti.^^^ results in differ- 
ent spin splittings [Fig.IHc)] of -25{^/i?- 2j/Re[,5Ke''']/i? 
and 2S^/R — 2j/Re[(5Ke'^]/i? for the conduction and va- 
lence bands, respectively. This explains the electron-hole 
asymmetry observed in the recent experimenti^. Here we 
remark that neither the off-diagonal {6k) nor the diago- 
nal ((5J^) term of the SOC alone can generate the electron- 
hole asymmetry since the two spin splittings can differ by 
sign at best, which actually implies the same magnitude 
of the spin splitting (see Fig. [1] for the sign convention) . 
Thus the interplay of the two types is crucial for the 
asymmetry. 



III. THEORY AND MODEL 

We calculate the 6k and (5^ analytically using degener- 
ate second-order perturbation theory and treating atomic 
SOC and geometric curvature as perturbation. For sim- 
plicity, we evaluate 6k and 6'^ in the limit k = 0. Al- 
though this limit is not strictly valid since k = does not 
generally satisfy the quantization condition on k^, one 
may still take this limit since the dependence of 6k and 
6'^ on k is weak. An electron at the K point is described 
by the total Hamiltonian -|- Hgo + He, where He 

describes the curvature effects and H^-^^^ describes the 
TT and a bands in the absence of both Hso and H^. The 
TT band eigenstates of iJ^'C^) are given by 
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where vp is the Fermi velocity and the momentum com- 
ponent kx along the circumference direction satisfies the 
quantization condition kx = (1/3-R)j^ for a {n,m) CNT 
with n — m = Sq + v{q'E'Z and v = ±1, 0) and 9 = 
tan~^[-\/3m/(27i + to)]. For a semiconducting {i> = ±1) 
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I* 



K,(0)v 

T(i) ' 



-i=(^e-«|V^^)±|^g))xT(i) (5) 



0. 



with the corresponding eigenvalues E^i^^^'^ = E^'^'' 

Here |^E'j^^^^■') with the upper (lower) sign amounts to the 
k = limit of the eigenstate at the the conduction band 
bottom (valence band top). XT(i) denotes the eigenspinor 
%B)) = ^E.=r.,„.e'''"-|PP is the orbital 

K,(0)\ 



of ay. IV-Afs)) 



projection of l^"-. 



^^•^ I into the sublattice AiH), \p^^ rep- 
resents the Pz orbital at the atomic position r, and the z 
axis is perpendicular to the CNT surface. 

When Hso and H^ are treated as weak perturbations, 
the first order contribution H^q to the effective SOC van- 
ishes since it causes the inter-band transition (Fig. [3]) 
to the (7 band^. The next leading order contribution to 
the effective SOC comes from the following second order 
perturbation Hamiltonian H^'^'^^—, 



FIG. 1: (Color online) Schematic diagram of the lowest con- 
duction (red, -E > 0) and highest valence (blue, _E < 0) band 
positions of a semiconducting CNT predicted by Ti^i^^^+Ti^oc 
for (a) (5k = ^k = 0, (b) 5k /O, 5k = 0, and (c) 5k /O, 5k /O. 
In (c), the conduction or valence band has larger spin split- 
ting depending on the sign of v. Arrows (green) show the 
spin direction along the CNT. The expressions for the energy 
level spacing are also provided. When they are negative, the 
positions of the two spin-split bands should be swapped. 



-pK 



EiO) - ijK,(0) 



H.C., 



(6) 



where the projection operator is defined by V^ = l — 
E.=T,il*"'^°^)(*^'^"^l- Another spin-dependent second 

order term H^oiV^ /{E^^^ -H^-^°^)]Hs^ is smaller than 
Eq. ^ (by two orders of magnitude for a CNT with 
i? ~ 2.5 nm), and thus ignored. Then the second order 
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energy shift E^(^l^^^ is given hy^ 



E 



K,(2) 



-E- 



K,(2) 



(7) 



where the upper (lower) sign appUes to the energy shift 
of the conduction band bottom (valence band top) and 
(V'5|i?'^'(2)|i/'^) = (V'i|i?'^'^^^|V'i) is used. Then by com- 
paring E^f^l^'^^ with Fig. [H one finds 

|=(^«|i?«^(^)|^i), | = (^^|7?^^(^)|^S). (8) 

Note that Sk and are related to pseudospin-flipping 
and pseudospin-conserving processes, respectively. 

To evaluate Eq. ([8]), one needs explicit expressions for 
Hso, He, and H^^^ . H^o is given by XsoY^r'^r ■ Sr^, 
where and Sr are respectively the atomic orbital and 
spin angular momentum of an electron at a carbon atom 
r. The tight-binding Hamiltonian of the Hs„ can be 
written^ as H,o = (Aso/2) E,.=r. (c^^ 



■ -st y 



'■A/ 1 

■ y\ X 



^) -f H.C., where 

<'r-i-(-)' '"r+(-)' ^'^'^ ^enotc the annihilation opera- 

tors for \pI)x+(-), bpx+(-), and \pl)x+(-)- Here 
denotes the eigenspinor of (+/— for outward/inward). 
For later convenience, we express X+(-) in term of X]{i) 
to obtain a expression for Hso, 



H.c. 



(9) 



For the curvature Hamiltonian He, we retain only the 
leading order term in the expansion in terms of a/ R. Up 
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FIG. 2: (Color online) Two-dimensional honeycomb lattice 
structure. x(y) is the coordinate around (along) the CNT 
with chiral vector na.i+ma.2 = (n, m) and chiral angle 9. 
ujj (j = 1, 2, 3), the length between y axis passing A atom and 
its parallel (red dashed) line is related with by ~u!j/{2R) 
[Eq. (Uni)]. The coordinates for the CNT is illustrated on the 
right. Here, x = ipR. 



to the first order in a/R, He reduces to H^"^ , 
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-I- V (r^^ — r-y^ 



(10) 



H.c. 



where is a lattice site in the sublattice A and its 
three nearest neighbor sites in the sublattice B are rep- 
resented by Bj (j = l,2,3) (Fig. 12). Here Sj, Xj, are 
proportional to a/R and denote the curvature-induced 
coupling strengths of s, Px, Py orbitals with a nearest 
neighbor pz orbital. Their precise expressions that can 
be determined purely from geometric considerations, are 
given by Sj = £,jSj, Xj = ijXj, and Yj = S^jYj with 
Ci w a/(2V3i?)sin6l, 6 ~ a/(2^/3i^) sin (7r/3 — 6), and 
6 ~ a/{2ViR) sin (tt/S + 0)(Fig. E]). Here 
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KpSinf 
K^cos 



6+^)' 



CpCos(J^ 

-v;^sinH-v;^~v;^cosH, 



pp 



Vpp sin- 



- y 
pp 



y-^ - I/O- 
sin(20)^V^, 



Sin 20--) 

3/ 2 



(11) 



Lastly, for the factor we use the Slater-Koster 

parametrizatior>i^ for nearest-neighbor hopping. In a 
band calculation, s, px, and py orbitals are used as basis. 

Combined effects of the three factors H^o, 'P^/{E^°'> - 
//K,(o)-j^ iJc are illustrated in Fig. [3] The real spin de- 
pendence arises solely from Hgo, which generates the 
factor CTj/^. For the pseudospin, the combined effect 
of Hso and H^ is to flip the pseudospin. When they 
are combined with the pseudospin conserving part of 
7?K/(-£;(o) _ ^K,(o))^ one obtains the pseudospin flip- 
ping process [Eq. ([8])] determining Sk- In addition, 
/{E^^'^—H^''--^^) contains the pseudospin flipping part, 
which is natural since states localized in one particular 
sublattice are not cigenstates of H^'^^K When the pseu- 
dospin flipping part of V^/iE^^^f — iJ^'(o)) is combined 
with Hso and He, one obtains the pseudospin conserving 
process [Eq. ([8])] determining (5^. 

The signs of SkG^^ and S'^/ cos 39 are negative. We 
find |((5{^/cos36')/(5k| = 4.5 for tight-binding parame- 
ters in ReffSi. Thus is of the same order as S}<r^, 
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FIG. 3: (Color online) Schematic diagram of the second order 
transition process generated by H^'^^^ [Eq. (|6}]. Pseudospin 
transitions (between the sublattices A and B) and interband 
transitions (between n and a bands) are illustrated. 



which is understandable since pseudospin flipping terms 
in £■(")- ffK,(o) ^^j^jj amplitudes V^p, Vfp) are compara- 
ble in magnitude to pseudospin conserving terms (with 
amplitudes i?^"^ — e^i-j,)). 



IV. BEHAVIOR IN A MAGNETIC FIELD 

Next we examine further implications of our result in 
view of the experimenl^, where the conduction band 
bottom and valence band top positions of semiconduct- 
ing CNTs {i' = ±1) are measured as a function of the 
magnetic field B parallel to the CNT axis. We find that 
the 9 dependence [Eq. ^] of S'^^ has interesting impli- 
cations. When cos 30 is sufficiently close to (close to 
armchair- type) , \S^\ is smaller than |(5Ke*^|. The predic- 
tion of our theory in this situation is shown in Figs.Ua) 
and (b) . Note that the spin splitting of both the conduc- 
tion and valence bands becomes smaller as the energy 
E increases. On the other hand, when cos 36* is suffi- 
ciently close to 1 (close to zigzag-type), \S'^\ is larger 
than |i5Ke'^|. In this situation [Figs. [D^c) and (d)], the 
energy dependence of either valence or conduction band 
is inverted; For i' = +l{—l), the spin splitting of the va- 
lence (conduction) band becomes larger as the the energy 
increases. 

Combined with the electron-prevailing [Figs.jDJa) and 
(c) for i/ = -|-l] vs. hole-prevailing [Figs.jH^b) and (d) for 
i/ = — 1] asymmetries in the zero- field splitting, one then 
finds that there exist four distinct patterns of E vs. B 
diagram, which is the second main result of this paper. 
Among these 4 patterns, only the pattern in Fig. [H^a) is 
observed in the experiments^, which measured two CNT 
samples. We propose further experiments to test the ex- 
istence of the other three patterns. 

Here we remark that although Eqs. H]), ([2]), ([3]) are 
demonstrated so far for semiconducting CNTs, they hold 
for metallic CNTs {v = 0) as well. For armchair CNTs 
with cos 30 = 0, 6^ becomes zero and the spin splitting is 



determined purely by Sk- For metallic but non-armchair 
CNTs, finding implications of Eq. ([1]) is somewhat tech- 
nical since the curvature-induced minigap appears near 
the Fermi level^^. Our calculation for (37, 34) (cos 36* «0) 
and (60, 0) (cos 36* = 1) CNTs including the minigap effect 
indicates that they show behaviors similar to Fig. IHb) 
and (d), respectively. Thus nominally metallic CNTs ex- 
hibit spin splitting patterns of iy = —l CNTs. 
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FIG. 4: (Color online) Calculated energy spectrum of the 
conduction band bottom (red, _E > 0) and valence band top 
(blue, E < 0) near K (solid lines) and K' (dashed lines) points 
in semiconducting CNTs with R « 2.5 nm as a function of 
magnetic field B parallel to the CNT axis. The chiral vectors 
for each CNT are (a) (38,34), (b) (39,34), (c) (61,0), and (d) 
(62,0), respectively. Arrows (green) show spin direction along 
the CNT axis and Ago denotes the zero-field splitting. Assum- 
ing ky = 0, the energy E including the SOC, the Aharonov- 
Bohm fiusi^ 4>ab = BtxE? , and the Zeeman coupling effects is, 
£; = ±ft«FA/(fc- + (l/^)(0AB/</)o))2-FS^7iT''"'+(^?/2)/^Br||B, 
with upper (lower) sign applying to the conduction (valence) 
band. (t)o = hc/\e\, =+!(-!) for XT(i)> vf ^ 
and g = 3^^. For estimation of E^^^^^'^'^\ we use tight- 
binding parameters in Ref.— ; Vss ~ — 4.76eV, V^p = 4.33eV, 
Vp'^ = 4.37eV, V^p = -2.neV, e« = -6.0eV, and ep = i^. 
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V. DISCUSSION AND SUMMARY 

Lastly we discuss briefly the effective SOC in a curved 
graphenei^. Unlike CNTs, there can be both convex- 
shaped and concave-shaped curvatures in a graphene. 
Wc first address the convex-shaped curvatures. When 
the local structure of a curved graphene has two princi- 
pal curvatures, l/i?i and l/i?2 with the corresponding 
binormal unit vectors rii and 112, each principal curva- 
ture l/Ri (z = l,2) generates the effective SOC, Eq. ([T|), 
with ffy replaced by cr • and R hy Ri. The correspond- 
ing Si and S'l values are given by Eqs. ([2|) and ([3|) with 9 
replaced by 9i, where 9i is the chiral angle with respect 
to n;. Thus the diagonal term of the effective SOC is 
again comparable in magnitude to the off-diagonal term. 
For the concave-shaped curvatures, we find that the two 
types of the SOC become —Si and —(5- with 9i, respec- 
tively. We expect that this result may be relevant for the 
estimation of the spin relaxation length in graphenes^ 
and may provide insights into unexplained experimen- 
tal data in graphcnc-bascd spintronic systems^!. We 
also remark that the effective SOC in a graphene may 
be spatially inhomogencous since the local curvature of 
the nanometcr-scalc corrugations^^ is not homogeneous, 
whose implications go beyond the scope of this paper. 

In summary, we have demonstrated that the interplay 



of the atomic SOC and the curvature generates two types 
of the effective SOC in a CNT, one of which was not 
recognized before. Combined effects of the two types of 
the SOC in CNTs explain recently observed electron-hole 
asymmetric spin splitting^ and generates four qualita- 
tively different types of energy level dependence on the 
parallel magnetic field. Our result may have interesting 
implications for graphcncs as well. 

Note added. -While we were preparing our manuscript, 
we became aware of a related paper—. However the ef- 
fective Hamiltonian [Eq. ([I])] for the SOC and the four 
distinct types of the magnetic field dependence (Fig. |4]) 
arc not reported in the work. 
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